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Abstract. Simultaneous eigenfunctions of two Askey- Wilson second order 
difference operators are constructed as formal matrix coefficients of the prin- 
cipal series representation of the modular double of the quantized universal 
enveloping algebra Wq{s[2{C)). These eigenfunctions are shown to be equal 
to Ruijsenaars' hypergeometric function under a proper parameter correspon- 
dence. 



1. Introduction 

The main goal of this article is to construct a solution of two commuting Askey- 
Wilson second order difference equations using representation theory of the modular 
double of the quantized universal enveloping algebra of 512(C). Furthermore we 
relate this solution to Ruijsenaars' hypergeometric function from jllj. 

By |S] there exist three inequivalent ^-structures on Uq, one associated to the real 
form su(2) of s[2(C), one associated to 5u(l,l), and one to s[2(]K). Koornwinder 
Noumi & Mimachi 9 , and Koelink 4 have shown that the Askey- Wilson 
polynomials arise as matrix coefficients of ^-unitary irreducible representations of 
Uq{B\x{2)). To prove these results they used the fact that the Askey- Wilson second 
order difference operator arises as the radial part of the quantum Casimir in lAq when 
calculated with respect to Koornwinder's twisted primitive elements. In jK] and 
[TF Koelink and Stokman constructed the trigonometric Askey- Wilson functions as 
matrix coefficients of ^-unitary irreducible representations of ZYg(su(l, 1)). In this 
paper we consider matrix coefficients of Z^g(sl2(R))-representations. 

An essential tool is the embedding of Z^g(s[2(M)) in Faddeev's |2 modular double 
of lAq. The modular double consists of two commuting copies of the quantized 
universal enveloping algebra of s[2(C) with deformation parameters q = 
and q = e^™^/^^ {wi,W2 S R>o), respectively. Kharchev et al. ^ made the crucial 
observation that the algebraic version nx of the principal series representation of 
Uq{sl2{M.)) on the space Ai of meromorphic functions on C can be extended to a 
representation of the modular double on the same space. In the same article they 
construct generalized Whittaker functions as matrix coefficients of tt\ . 

We construct joint eigenvectors to the action under tt\ of two commuting twisted 
primitive elements (one for each copy of the quantized universal enveloping algebra 
of s[2(M) inside the modular double) in terms of Ruijsenaars' 10 hyperbolic gamma 
function. The action of the two commuting quantum Casimir elements in the 
modular double shows that the corresponding matrix coefficients, for which we have 
an explicit integral representation, satisfy Askey- Wilson second order difference 
equations in step directions iwi and iw2 ■ By construction these matrix coefficients 
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are invariant under interchanging of wi and W2 ■ We show duality of this solution S 
in its spectral variable A and its geometric variable. Consequently it satisfies another 
two Askey- Wilson second order difference equations in its spectral variable. 

In a series ^Sj of papers, Ruijsenaars considered a solution R of the 

same Askey- Wilson difference equations. These equations arose in his study of 
relativistic quantum integrable systems. Ruijsenaars defined the hypergeometric 
function i? as a Barnes' type integral with integrand expressed in terms of the hy- 
perbolic gamma function. Subsequently he established for R duality, D4-symmetry 
in the parameters, asymptotic behaviour and the reduction to Askey- Wilson poly- 
nomials. We use these properties to show equality of R to S, which is not apparent 
from their explicit integral representations. 

The structure of this paper is as follows. In Sections |21 and 13 we recall some 
properties of the hyperbolic gamma function and of Ruijsenaars' hypergeometric 
function R, respectively. In Section 0] we define the modular double of Uq and its 
principal series representation on meromorphic functions. In Section Owe consider 
the corresponding eigenvalue problem of two commuting twisted primitive elements. 
Using the matrix coefficients of the principal series representation we construct a 
solution S to the Askey- Wilson difference equations in Sectional and we establish 
the duality of S. In Section [7| we show by a direct calculation that S reduces to 
the Askey- Wilson polynomials for certain discrete values of the spectral parameter. 
Finally, in Section Owe show that S equals Ruijsenaars' hypergeometric function 
R. 

1.1. Notational conventions. If ± appears inside the argument of functions we 
mean a product, e.g. 

(1.1) f{z±a)^ f{z + a)fiz-a). 

Otherwise it means that all sign combinations are possible. 

Whenever we use a square root, we always mean the branch which has a cut 
along the negative real line and maps the positive real line to itself. 

1.2. Acknowledgements. I would like to thank my supervisor Jasper Stokman 
for many illuminating conversations and his help in writing this article. 

2. The hyperbolic gamma function 

Both Ruijsenaars' and our solution to the Askey- Wilson second order difference 
equations are expressible in terms of the hyperbolic gamma function, which was 
introduced in ^Ul- Let us therefore recall some basic properties of this function, 
see and the appendices of [U for more details. 

Let us first define for wi,W2 G C+ = {z E C \ 5R(z) > 0}, 

sin(2j/z) z \ dy 

2sinh(wij/) sinh(w22;) wiW2y J y 

Notice that the integrand has no pole at 0. To ensure convergence of the integral at 
infinity however, we must impose the condition |3(2)| < 5R(w), where w is defined 

by 

Wi + W2 
W ~ . 
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The hyperboUc gamma function G{z) = G{wi,W2', z) for < is now 

defined by 

(2.2) G{wi,w2;z) ^€'^^'^^•'"'•'1 

The hyperboHc gamma function G owes its name to the fact that it satisfies the 
difference equations 

G{z + iwi/2) = 2cosh(7rz/w2)G(z — i'Wi/2), 

"^^■^^ G{z + iw2/2) ^ 2cosh(7rz/u;i)G(z - ^2/2). 

In these equations we suppress the wi and W2 dependence of G, which we continue 
to do whenever this does not cause confusion. These two difference equations allow 
for an analytic continuation of G to a meromorphic function on C. The hyperbolic 
gamma function can also be expressed in terms of Barnes' double gamma function, 
or Kurokawa's double sine function. Details can be found in ^] Appendix A]. 

Let us first note a few symmetries of the hyperbolic gamma function, which are 
all obvious from H2.1() : 

(2.4) G{wi,W2]z) ^ G{w2,wi]z), 

(2.5) G{wi,W2\z) ^ G{wi,W2\-z)~'^, 

(2.6) G(wi,W2\z) = G{wi,w2,-z), 

(2.7) G{p.Wl,^JLW2\^J.z)^G(wl,W2\z) (^eR>o). 

The pole and zero locations of G are easily derived from the difference equations 
(12.31) . since G has no poles or zeros in the strip z e R x i(— 5R(w), 'R[w)) in view of 
(12.2(1 . The zeros of G are contained in the set 

(2.8) A+ — iw + iwiZ>o + iw2'^>a 

and the poles in — A+. The pole at z — —iw is simple, and its residue equals 

(2.9) T^V^i^- 

ZTT 

If w\lw2 is irrational all other poles are also simple and their residues can be 
calculated from (|2.9() and the difference equations (|2.3(l . see [T^ Proposition III. 3]. 

For later purposes it is convenient to call an infinite sequence of points in C 
increasing (respectively decreasing) if it is contained in a set of the form a + A+ 
(respectively a — A+) for some a S C. In this terminology, G has one increasing 
zero-sequence and one decreasing pole-sequence. 

We also need an estimate for G(z) as 5ft(z) 00 and 3(z) stays bounded. In 
fact we only need it for the quotient of two hyperbolic gamma functions, which is 
easily derived from the estimate of the hyperbolic gamma itself as described in jlUI 
Proposition III. 4] and |12[ (3.3)]. For a,b £ C and wi,W2 G (0,oo) the resulting 
estimate reads 

G{z~a) / TT 2 .2 



(2.10) -( = exp (22(6 -a) + a'-b' + f{z)) 

G[z — 0) \2iwiW2 

where f{z) satisfies for di{z) > max(wi,W2) + max(3fj(a), 5R(6)), 

(2.11) |/(z)| < G(«;i,w;2,3(^),a,6)e-"^(")/'""'^('"i'"'^), 
with C depending continuously on (0, 00)^ x M x C^. 
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We also use the description of G as a quotient 



(2.12) G(z) 



where E is an entire function with zeros at A_|_ which are all simple if 'Wi/w2 is 
irrational. For a precise definition of E, see ^ Appendix A]. 

We will occasionally meet functions defined by an integral of the form 

for wi,W2 > and for parameters Uj and dj satisfying |3(uj)|, |3(dj)| < w and 
^(Sj=i i'^j^^j)) > 0- These conditions ensure that the integral is well defined (the 
contour meets no poles and it decreases exponentially at ±oo). In '11' Appendix 
B] it is shown that 

n 

M{u, d) Y[ E{-iw + Uj - dk) 

3,k=l 

has a unique analytic extension to the set {{wi,'W2,u,d) e x C^" | 3(X]("j ~ 
dj)/wiW2) > 0}. Hence M{u,d) is a meromorphic function which can only have 
poles when some E[—iw + Uj — dk) is zero. 



3. RUIJSENAARS' HYPERGEOMETRIC FUNCTION 

Ruijsenaars JI] introduced a generalization R of the hypergeometric function as 
a Barnes' type integral. We recall several properties of R from and ^21 which 
we will need to relate R to the formal matrix coefficients we are going to define in 
subsequent sections. 

We define Ruijsenaars' hypergeometric function in terms of a parameter set 7^ 
(^ = 0, 1, 2, 3), which is related to Ruijsenaars' original c-parameters by (l-H)]- 
Dual parameters 7^ are defined as 



(3.1) 
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We denote the set of parameters (70, 71, 72, 73) by 7 and the set of dual parameters 
by 7. Note that taking dual parameters is an involution, 7 = 7. 

Ruijsenaars' hypergeometric function R is now defined for generic parameters 

wi,w2 e C+ , 7 e by 

(3.2) Rir,x,X) = 

1 P G{z ±x + ijo)G{z ± A + 170) T-r G(i7o + ijj + iw) 
^/wlW2 Jc G{±x + i7o)G(±A + «7o)G(z + iw) j-J-^ G{z + 470 + i"fj + iw) 

Note that we use the convention Hl.l(l in this expression. The integral is taken over 
a contour C, which is a deformation of M seperating the poles of the numerator 
from the zeros of the denominator (equivalently, C separates the increasing pole 
sequences of the integrand from the downward pole sequences). R has an analytic 
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extension to a meromorphic function on (wi, W2, 7, 2:, A) G x C^, with possible 
poles for fixed values of wi, W2, and 7 at 

(3.3) xe±(K+-ii,), Ae±(A+-i7,) (j=0,l,2,3). 

Recall that A+ is defined by (|2.8|) . 

We now look at the Askey- Wilson second order difference equations which R 
satisfies. The equations are obtained from (lli Theorem 3.1] by not only replacing 
the c- variables by 7, but also multiplying the equations by a constant. These 
descriptions of the Askey- Wilson difference equations are more convenient for the 
representation theoretic approach we consider in the following sections. 

Let us define the function A by 

(3.4) 

g7ri«)i/i02+2iri7o/i"2 JL / ly^^ 

A{wi,W2,i;x) = J- — r , , M cosh — {x + — + ilj) 

sinh(^)sinh(Hl(^) -L^o V^«2 2 



f 
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The Askey- Wilson second order difference operator is defined by 

(3.5) C^f{x) = A{wi,W2,j;x){f{x + iwi) - f{x)) 

+ A{wi,W2,'^;-x){f{x-iwi)-f{x)). 

Here the superscript x is added to emphasize that the operator acts on the x- 
variable (in a moment we will also consider the operator £ acting on the spectral 
variable A). We write for the Askey- Wilson operator H3.5|) with wi and W2 
interchanged. 

Ruijsenaars' hypergeometric function R is an eigenfunction of four Askey- Wilson 
second order difference operators with eigenvalues expressible in terms of 

(3.6) v{wi,w2,j;X) = 

_ 2e'^™i/'"^+2'^''^»/"'=^(cosh(27rA/w2) + cosh{7riwi/w2 + 2Tri%/w2)). 

Specifically, R satisfies the difference equations 

£^i?(7; X, A) = v{wi,w2,j; X)R{-f; x, A), 

C%R{-i] X, A) = v{w2,wi,j; X)R{^; x, A), 
(3.7) 

£^i?(7; X, A) = v{wi,W2,'j; x)R{'^; x, A), 
C^R{j; X, A) = w(w2, wi, 7; 2^)^(7; x, A). 

Actually the last three of these equations follow from the first by various symmetries 
of R. The second difference equation can be derived from the first (and the fourth 
from the third) by using the fact that R is invariant under the exchange of wi and 

W2, 

R{wi,W2,^;x, X) = R{w2:Wi,j;x,X). 

This symmetry can be directly seen from the definition 1)3. 2|l of R and the corre- 
sponding symmetry 1)2. 4|1 of the hyperbolic gamma function. The third difference 
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equation can be obtained from the first by using the duaUty of R under the exchange 
of X and A, 

(3.8) R{r,x,\) = R{fA:x). 

This duahty is also a direct consequence of the definition of R using the fact that 
70 + 7j = 70 + 7j for j = 1, 2, 3. 

There are more symmetries of R directly visible from the definition. Since the 
hyperbolic gamma function is scale invariant it follows that R is scale invariant as 
well, 

R{uwi, VW2, v^; vx, vX) — R{wi,W2, 7; x, A) 

for u e (0, 00), where denotes the scaled parameter set (1^70, 1^71, 1^72, ^^73)- Fur- 
thermore it is immediately clear that R is symmetric under permutations of 71, 72, 
and 73 . This symmetry can be extended to a £'4-symmetry in the four parameters 
7 (where the Weyl group of type -D4 acts on the parameters by permutations and 
an even number of sign flips). To formulate this result we need the c- function 

1 ^ 
G{2y + iw) j-J-^ 

and the normalization constant 

3 

(3.9) N{J) = l[G{iJo + ^7J+^w). 

The il'4-symmetry ^1 Theorem 1.1] of R then reads 

R{'i;x,X) R{w{'y);x,\) 



(3.10) 



c(7; x)c{f, X)N{j) c{w{-f);x)c{w{j); X)N{w{-f)) 



for all elements w of the Weyl group of type D4. Notice that both the c- function 
and N are invariant under the action of the S'3-subgroup which permutes 71, 72, 
and 73. 

Finally we recall the limit behaviour of R, cf. |12l Theorem 1.2]. Set a = 
2tt/wiW2- For wi,W2 S K>o, 7 G M"*, and wi ^ W2 there exists an open neighbour- 
hood U C C of M, such that the asymptotics of R for fixed A € J7 are given by 



(3.11) i?(7; X, A) ^ c'(e"(ia(^)l-T'o~"')l3^(^)l) 

for ^{x) ±00, uniformly for 3(a;) in compacta. In fact, Ruijsenaars gives a 
precise expression for the leading term of R as '3i{x) — > ±00 when A S M. These 
results easily extend to A in some open neighbourhood U of R. 



4. The modular double of Uq{5l2{C)) 

In this section we consider a slightly extended version of Faddeev's [2] modular 
double of Uq{sl2{C)) and define an algebraic version of its principal series repre- 
sentation on the space M of meromorphic functions on C. We define an inner 
product on some suitable subspace of A4, which is compatible to the ^-structure 
on Uq{sl2{C)) associated to the real form s[2(K) of s[2(C), cf. jH). 
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Throughout Sections 01 [3 and El we assume that wi and W2 are positive real 
numbers such that neither wi/w2 nor 1x12 /wi is an integer, unless specifically stated 
otherwise. We define 

g = exp(7riK;i/zi;2), g = exp(7riw2/wi), 

which both lie on the unit circle (but they are not ±1). For complex numbers j3 we 
define 

Definition 4.1. The quantized universal enveloping algebra lAq = ^^^(5(2(0)) 0/ 
s[2(C) is the unital associative algebra over C generated by K^^ , E, and F , subject 
to the relations 

KK~^ ^ K-^K = 1, 
KE = q^EK, 
KF = q^^FK, 

EF^FE=^-^-' 



q-q ^ 



If Wi/w2 is irrational, then the center of Uq is generated by the quantum Casimir 
element 51, defined as 

n = qK + q-^R-^ + {q - q-^fFE. 

By simply replacing qhy q (or interchanging wi and W2) we obtain the quantum 
universal enveloping algebra Zig. The generators oiUq are denoted by , E, and 
F. The following concept of modular double was introduced by Faddeev |21. 

Definition 4.2. The modular double Q islAq^Uq endowed with its standard tensor 
product algebra structure. 

For elements X eUq (respectively X G Uq) we also write X (respectively X) for 
its image under the natural embedding of Uq (respectively Uq) in Q. In particular, 
XX = XX in Q for elements X e Uq and X eUq. 

We now define an extension of the modular double by formally adjoining com- 
plex powers of K and K to Q. Let A = ®xi^c'^^ ^'^ group algebra of the 
additive group C = (C, ©), where © is the translated addition x(By = x+y+iw (this 
translation in addition will make formules simpler later on). The unit of C is —iw. 

Lemma 4.3. There exists a unique left A-action by algebra automorphisms on the 
modular double Q satisfying 

X- K^^ ^ x-K^^^K^^, 

X- E = -qe^'"'''"^E, x-i^ = -ge^^^/^ii^, 

x-F^ -qe-2^^/'"^F, x-F = -qe-2'^^/"'iF. 

Proof. Observe that e.g. the action of a; on £' can be rewritten as 

X ■ E = ^27r{x+iw)/w2-g_ 

The lemma now follows by direct calculations. □ 
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Definition 4.4. The extended modular double D = Qyi A is the crossed product of 
the modular double Q and the algebra A under its action on Q as defined in Lemma 

113 

Hence T) is the vector space A endowed with the unique algebra structure 
such that the natural embeddings of Q and AinV are algebra morphisms and such 
that 

xQ={x- Q)x, Wx e C, VQ e Q, 

where we identified x (respectively Q) with their images under the natural embed- 
dings of A (respectively Q) in V. 

Now we define representations tt\ of the extended modular double Q on the space 
M of meromorphic functions on C depending on a complex representation label A, 
cf. in]. These representations may be viewed as algebraic versions of the principal 
series representations oiUq{sl2{M.)). We define these representations in terms of the 
operators Ty and Sy on A^, which act by 

Tyf{z)^f{z + y), Syf{z)^e^-^^'yf{z) (y G C). 

Lemma 4.5. For A G C the assignments 

,1/2 



q-q 

uniquely define a representation Tr\ of V on M. 

Observe that the action of the generators of Uq are obtained from the action of 
the generators of Uq by interchanging wi and W2 ■ 

Proof. The defining relations of T) are easily checked using T^Ty = T^+y ~ TyT^, 
SxS-x = 1, and the equation 

-^x'^y ^ '^y-^x- 

□ 

Remark 4.6. Denote v = {wi —W2)/2, then Tixiiv) — nx{K) and nx{—iv) — n\{K). 
The extension of the modular double Q by ^ and the extension of the representation 
tt\\q to tta thus have the effect of introducing non- integral powers of Tiu^^ and T^^^ 
in the image of tt\. The introduction of this extension is not an essential part of the 
analysis later on and is only included for simplification. Using only integral powers 
of K and K we can simulate the action of x for x in some dense subset of M, cf. |S1 
Proposition 1.6]. 

A simple calculation shows that 7T\{n) acts as a scalar, 

(4.1) 7rA(r!)/ = -2cosh(2^AM)/, f e M. 



RUIJSENAARS' HYPERGEOMETRIC FUNCTION AND THE MODULAR DOUBLE 9 



Since tt\ is an algebraic version of the principal series representation with represen- 
tation label A G C, this is as expected. 

Definition 4.7. We say that f Cz A4 has exponential growth with growth rate 
e e K i/ there exists a compact set Kf € R such that all poles of f are contained in 
Kf X iR — {x + iy \ X £ Kf,y £ R} and if \ f{x + iy)\ = 0{exp{e\x\)) for x ±00, 
uniformly for y in compacta of R. 

On the space of mcromorphic functions which have negative exponential growth 
and which have no poles on R, we define a sesquilinear form by 

(4.2) / f{z)W)dz. 

Observe that this expression is already well defined under the milder asymptotic 
condition that the sum of the two exponential growths of / and g is negative. Note 
furthermore that H4.2|l can be rewritten as 

(4.3) {f,g) = f f{z)-g{z)dz, 

Jr 

where g(z) := g{z) now is analytic at z e R. 

Following ^ we define an antilinear anti-algebra involution * on the extended 
modular double T) by 

(4.4) K* ^K, E* ^ -E, F* = -F, K* =K, E* = ~E, F* = -F, x* = ^x. 

If we restrict this involution to lAq (respectively lAq) we obtain the ^-structure on 
Uq (respectively Uq) corresponding to the noncompact real form s[2(M) of s[2(C), 
cf. 0. 

The following lemma relates the sesquilinear form H4.2|l to the ^-structure H4.4|l 
on T). 

Lemma 4.8. Let A G C and f,g £ A4. If the poles of f and g are outside the strip 
M X i[—wi, wi] and if the sum of the exponential growth rates of f and g is smaller 
than ~2tt/w2, then 

{TT^{X)f,g)^{f,TT-^{X*)9) 

for X £ Uq.i := spanc{l,E, F, K, K-\ FK, ER-^}. 

Proof. In view of H4.3() the proof follows by a change of variables and some contour 
shifting using Cauchy's theorem. □ 

A similar lemma holds for the dual algebra Uq. 

5. Twisted primitive elements and matrix coefficients 

Koornwinder 7 introduced twisted primitive elements to obtain the Askey- 
Wilson polynomials as matrix coefhcients of finite dimensional Z-/q(s 12(C) )-repre- 
sentations. We recall the definition of twisted primitive elements and show that 
they act as first order difference operators under the representations tta. We con- 
struct eigenvectors to these operators in terms of the hyperbolic gamma function 
and we consider the corresponding formal matrix coefficients ip of tt^. In subsequent 
sections we relate tp to Ruijsenaars' hypergeometric function. 

Let p £ C and set 

Up = g2p/«.i _^ q-2p/un ^ 2cos{2ttp/w2). 
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The twisted primitive element Yp ^ Uq C V is defined as 

(5.1) Y„ = iq-^^^E + iq-^'^FK "-^{K - 1). 

q-q 1 

Analogously we define the twisted primitive clement Yp G hiq by interchanging 
wi and W2, viz. 

Yp - iq-^'^E + iq-^/^Fk - 1), 

where Vp = 2cos(27rp/wi). 
Denoting 

I . Vp-Vr , , , Vp-Vr 

tJ'T[P) = ^J-t(P) ^ 

q-q 1 q-q i 

we now have the following lemma. 

Lemma 5.1. The function 

ttX ( \_ G{z + \/2-3iw/2±iT) 
"'"^^^ " Giz-X/2-iw/2±ip) 

satisfies 

7rx{Yp)H^^^ = tir{p)H^,p, 
7T^{Yp)H^^p^Pr{p)H^,p. 

Proof. Since p is invariant under the exchange of wi and W2, it is sufficient to 
prove only the first eigenvalue equation. A calculation shows that TTx{Yp)f = p.T{p)f 
is equivalent to the first order difference equation 

cosh(^(z + A/2 - ■iiw/2 ± ir)) 

(5.3) f{z + ^w^/2)^ ^ L^_±^-2L fi^^ ^ ^^^/2). 

cosh(^(z — A/2 — iwl2 ± ip)) 

(The exact calculation can be found in Appendix^) Using the difference equation 
(12.31) for the hyperbolic gamma function it immediately follows that -ff^p satisfies 
the difference equation H5.3|l . □ 

Remark 5.2. For any one of the two equations H5.2|l there are infinitely many so- 
lutions (we can e.g. multiply a solution to the first equation by any iwi-periodic 
function). The crucial step in finding common solutions to both difference equa- 
tions is to rewrite the first difference equation in the specific form (|5.3() . Indeed, 
the resulting solution ^ in terms of hyperbolic gamma functions is invariant un- 
der interchanging wi and W2, hence it automatically satisfies the second difference 
equation. This is the main difference between our analysis and the one in |15j . 

Now let us consider the adjoint Y* , which is 

Y; = iq'/^E + iq-^/^FK + -^^{K - 1). 

Since /z,j((t)* = — /io((T), we are interested in solutions to the equation TT^(Y*)f = 
—fj,fj{a)f and the corresponding equation T:j^{Y*)f = —jlfj{a)f for the second com- 
ponent of the modular double. 
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Lemma 5.3. The function 

, Gjz+'X- iw/2 ± iv) 

"^"^^^ Giz-\ + iw/2±ia) 

satisfies 

Proof. The proof is similar to the proof of the previous lemma. □ 

We will need a few results on the analytic properties of the two functions p 
and F!: 



- V,(T' 



Lemma 5.4. The possible pole locations of p and ^ are at 

-A/2 ± ir + - A+, A/2±ip + iw + A+ 

and 

-A/2±w-A+, A/2±iCT + A+ 
respectively. Furthermore, p and F^ ^ have exponential growth with growth rates 
7r(29(A) — 2w) / W1W2 and 7r(— 25(A) — 2w)/wiW2, respectively. 

Proof. The proof follows directly from the zero/pole locations and asymptotics of 
the hyperbolic gamma function (see Section |2J). □ 

Define 

(5.4) C = max(|5R(p)|, |5i(a)|, |5R(t)|, |5R(z;)|) 
and 

(5.5) C = W2-e-|3(A/2)|. 

We assume that the parameters p, cr, u, r and the variable A are such that C > 0. 
For 13(2^)1 < C define 

(5.6) ijip,a,T,v;x,X) - {TTx{x)H^^p, F^J, 

which is well defined since the exponential growth of the integrand equals —2tt{wi + 
W2)/wiW2 < and the pole sequences of iT\{x)H^p and all stay away from the 
real line due to the condition |3(a;)| < C- Note that the increasing pole sequences 
of TT\{x)H^p and F^^^ are all located above the real line and the decreasing pole 

sequences are all located below the real line due to the shifted addition in C. 
Observe furthermore that the matrix coefficient tp is invariant under the exchange 
of wi and W2, cf. Remark 15.21 The function i/j will be related to Ruijsenaars' 
hypergeometric function R in Section |S1 

Using (|2.5|) . H2.6I) . and (|4.3|) we can write as 
(5.7) 

f G(z + x + \/2- iw/2 ± iT)G(z - A/2 - iw/2 ± ia) , 
G(z + X — A/2 + iw/2 ± ip)G{z + A/2 + iw/2 ± iv) 

which is of the form (|2.13() . It follows from the discussion at the end of Section [21 
that 

^'(7; X, A) = E{x ± ijQ)E{x ± iji)E{-x ± ij2)E{-x ± 173) 
^^•^^ X E{X ± t%)E{X ± i^i)E{~X ± i^2)E{~X ± ^73)1^(7; X, A) 
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has an entire extension to 

(5.9) O = {(wi,W2,p,cr,r,u,x,A) e C+ x C^}. 

Hence ^ is meromorphic on the same domain O. 

6. The Askey- Wilson difference equations 

We show that the formal matrix coefficient 4' (see Section satisfies a second 
order difference equation with step size iwi using a radial part calculation of the 
Casimir f2 with respect to twisted primitive elements. As a consequence a renor- 
malization S H6.10|l of V' satisfies an Askey- Wilson second order difference equation. 
Since S*, like '0, is invariant under exchanging wi and W2, we obtain a second dif- 
ference equation with step size iw2 ■ We furthermore show that S satisfies a duality 
in the geometric and spectral variables, and we derive various obvious symmetries 
of S". 

Let us start by establishing a correspondence between the set of parameters p, 
cr, r, and v and Ruijsenaars' parameter set 7 by 

(6.1) 7o = -p-|-o-, 7i=p + cr, 72 = -T-t;, 73 = r - u. 

Observe that 70 (see (jH.lll ') becomes 

70 = ^(70 + 71 + 72 + 73) = o- - 

We will also use the abbreviation 7 for the parameters (p, a, r, v). In particular we 
write "0(7; x, A) for (|5.6|l . Later we show that 1)6. l|l is the parameter correspondence 
which relates ip to Ruijsenaars' hypergeometric function. 

Now we perform a radial part calculation of the Casimir element 57 with respect 
to the twisted primitive elements (see (15.1(1 '). The result is stated in terms of the 
function A, see (13.4(1 . 

Lemma 6.1. We have 

xVl — xn{x) mod xUq,2{Yp — f^rip)) + {Ya — Pv{<^j)xUq^2, 

where Uq^2 ■= spanf^{l, K^^ , F} and fl{x) is defined as the element 
n{x) = B{x)K + C{x) + D{x)K-^ 

with coefficients 

B{x) = g"^74(7o,-7o,7i,-7i;2^), 

C{x) = (7-i+27o/»i L ^(^. _ A[r, -x) + l + g2-47o/-i 
D{x) 9^^(72,-72,73,-73; -a;). 

Proof. The proof involves a radial part calculation similar to the one performed in 
[T31 Proposition 3.3]. In fact we can use the calculation in [T3] using an embedding 
(j) of the extended quantum universal enveloping algebra Uq into the one in [15j , 
given by 

V ^Wl J 

(t){E) = -iKX+, cj){F) = iX-R-^. 
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A direct calculation gives (jj^n) = {q-q-^)'^n + 2 and 0(yp) = F2p/«,i (on the right 
hand side we use the f2 and Y from |E1 > which have a slightly different definition) . 
Note that in [T5| the radial part is calculated modulo a larger vector space. However, 
it is easily verified that the present smaller space suffices for the proof. □ 

Using this radial part calculation we can prove that ip H5.6|l satisfies a gauge 
transformed Askey- Wilson second order difference equation. 

Lemma 6.2. The function ijj{x) — ?A(7;x, A) satisfies the difference equation 

(6.2) -2 cosh(27rA/w2)V'(a;) = B{x)^p{x + iwi) + C{x)tp{x) + D{x)^{x - iwi), 

and a similar equation with wi and W2 interchanqed. These equations hold as 
identities between meromorphic functions on the domain O {see H5.9|l V 

Proof. Observe that by the symmetry of ijj in wi and W2 we only have to prove the 
difference equation (|6.2|1 . 

We first prove the lemma under restricted parameter conditions, which allow 
us to use the expression H5.6|l of -0 as a matrix coefficient of the D-representation 
tta- Using analytic continuation we can subsequently remove these parameter con- 
straints, cf. the discussion at the end of Section [S] 

Let us assume that wi,W2 > and that 

(6.3) W2 > 7wi +4e + 2|3(A)| +4|3?(a;)| 

holds. Then |3(a;)| < C, so V' is defined by H5.6|l (recall that ^ and ( are defined by 
(|5.4(l and l|5.5|l . respectively). By 14.f(l . 

(6.4) -2cosh(2^A/u;2)0(a;) = {nxixn)H^^p, F^^,) 
holds. By Lemma Ffi. f I there exist X, Z E Uq^2 such that 

(6.5) xn = xn{x) + xX{Yp - ^ir{p)) + (y^ - ^i^{a))xZ. 
Since T^xiYp — ^ir{p))H^ p = 0, we have 

(6.6) {iTx{xX{Yp - ^ir{p)))Kp, F^J = 0. 

The exponential growth of tt\{xZ)H^ p is at most the exponential growth of ^ 
plus 27r/w2 (due to the possible occurrence of an Siw2 factor in ■K\{Zy). The sum of 
the exponential growths of n\{xZ)H^p and F^^^ is at most — 27r/wi, hence strictly 
smaller than —2tt/w2, since the restrictions on the parameters imply that W2 > wi. 
Moreover the condition (|6.3|l implies that neither tt\{xZ)H^ p nor ^ has any 
poles in the strip M x i[— wi, wi]. Using Lemma lOl and the fact that Y^, G Uq^i, we 
thus obtain 

(6.7) {ttx{{Y„ - ti^{a))xZ)Hr,p,F^^„) = {■K^{xZ)Hr,p,TT-,,{Y: + pL^{a))F„,,) = 0. 
Combining H6.4|l . (|6.5|l . (|6.6|) . and (|6.7|l now yields 

(6.8) -2cosh(2^AM)^(a;) = {^ximKp^F^^^) = {7r^{xn{x))H^^p, F^^^). 

Furthermore, by Lemma |6.fl fremember that x + iwi and xK act in the same way 
under Tr\) we have 

(6.9) {Trx{xn{x))Hr,p., F^,^) = B{x)iIj{x + iwi) + C{x)^p{x) + D{x)iIj{x - iwi). 

The lemma for the restricted parameter conditions follows now directly from H6.8(l 
and jnini- □ 
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Using the function 

, , G(x + i72)G(a; + 173) 
G{x - ijo)G[x - 171) 

we can define a renormalization S oi ip as 

^'■''^ ^(^^^'')=V^A(7;x)A(7;A)- 

The function N H3.9|) is a convenient normaUzation factor when matching S to R 
in Section |S1 

Lemma 6.3. 5'(7; x, A) is meromorphic on O with possible poles at 

X = ±{iy-i^k), x^±{v-ijk), i^Q + i^n ^ -V - iw 
for V e K+, k ^Q, 1, 2, 3, and 1 = 1,2, 3. 
Proof. Using (|5.8|l and (|2.12() we can express S as 

Sirx X) = ^ii;x,X)N{j) 

' ' llLoE{±x + ijk)E{±X + t%) 

From this expression we can easily read off that the possible pole hyperplanes are 
as stated in the lemma (they have to be either poles of N{j) or zeros of one of the 
i?-functions in the denominator). □ 

Theorem 6.4. The function 5(7; x, A) is a simultaneous eigenfunction of the two 
Askey-Wilson type second order difference operators CH^ and {see (|3.5|l ) with 
eigenvalues v{X;wi,W2,j) and v{X;w2,wi,j) respectively, where v is defined by 

Proof. Note that A satisfies the first order difference equation 

cosh( — (x + 172)) cosh( — (a; + ^73)) 

A(x + iwi/2) = , ) rriFl ^A(x - zu-i/2). 

^ ' cosh(4(a:-*7o))cosh(^(a;-*7i)) ^ ' 

The desired eigenvalue equation H3.5|l for now follows immediately from Lemma 

To prove the result for the operator we note that S is symmetric in wi and W2, 

while interchanging wi and W2 transforms C to C We could also prove the second 
difference equation by repeating the argument for the first difference equation using 
the component Uq of the modular double. □ 

We continue the analysis of the eigenfunction S by proving its duality in the 
geometric variable x and the spectral variable A, similar to the duality H3.8|) for 
Ruijsenaars' hypergeometric function R. The duality transformation 7 — > 7 of the 
parameters (see (|3.1|l ^ is equivalent to interchanging p and v under the parameter 
correspondence Hfi.f (I : (p, a, t, v) — > {v, a, t, p). 

Theorem 6.5 (Duality). We have 

S{j;x,X) = S'(7; A,a:;) 

as meromorphic functions on O. 
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Proof. Assume that wi,W2 > and w/2 > S, + |3(a;)| + |3(A)|, where is as 
in (|5.4|l . Note that these restrictions on the parameters are invariant under the 
exchange {x,"/) <-> (A, 7). Then we can use the integral representation (|5.7|l for 
both ipil'j X, A) and ■0(7! A, x) to compute 

f G{z + x + \/2~-iw/2± iT)G{z - A/2 - w/2 ± ia) 
X, A) - + a; ^ A/2 + iw/2 ± ip)G{z + A/2 + iw/2 ± iv) 

G{z + x/2 + A - iw/2 ± iT)G{z - x/2 - iw/2 ± ia) ^ 
G{z + x/2 + iw/2 ± 2^)0(2 - 2;/2 + A + iw/2 ± w) ^ 

where we used the change of integration variable z z + (A — a;)/2 and a contour 
shift in the second equality. This contour shift is allowed since the integrand con- 
verges to zero exponentially at ±c», and the conditions on the parameters ensure 
that there are no poles picked up by shifting the contour back to R. 

Since ^ (see (I5.8|l ') is entire on O, it follows that ■!/'(7;a^, A) — -0(7; A, a;) holds 
as identity between meromorphic functions on O. The desired duality for S now 
follows from N{'y) = N{j) and 7 = 7. □ 

Corollary 6.6. The function 5(7; x, A) is a simultaneous eigenfunction of the 
Askey-Wilson second order difference operators C^, C^, C^, and with eigenval- 
ues v{X;wi,W2,j), v{X;w2,Wi,j), v{x;wi,W2,^), and v{x;w2,Wi,^) respectively. 

Proof. The fact that S is an eigenfunction of £^ and £^ was proved in Theorem 
16.41 The proof for the other two difference operators follows from this fact and 
duality fTheorem l5.5|l . □ 

It is immediately clear from the integral representation (|5.7|) that ip is invariant 
under sign flips of the parameters p, a, r, and v. This leads to the following 
symmetries for S. 

Lemma 6.7. Let Wn be the Weyl group of type Z?„, which acts on n-tuples by 
permutations and even numbers of sign changes. Let V = W2 x W2 C W4 be the 
Weyl group of type D2 x D2, where the first (respectively second) component acts on 
the parameters (70,71) [respectively (72,73)) of the four-tuple (70,71,72,73)- For 
an element v G V we have 

S{'i;x,X) ^ S{v{-f);x,X) 

c(7; a;)c(7; A)iV(7) c{v{-f);x)c{v(j);X)N{v{-f)) 
as meromorphic functions on O. 

Proof. Note that the action of ^ ~ on the parameters {p,a,T,v) is by sign 
flips of p, a, T, and v. Under the conditions C > s.nd |3(a;)| < C it follows from 
the integral representation (|5.7|l of that is invariant under the action of V on 
7 (note that the parameter restrictions are ^-invariant). 
Observe that 

G(x±i72)G(a;±i73) 
<r,x)A{j;x)- G(2^^) 

is also invariant. Since the action of V commutes with taking dual parameters 
(which is obvious in the parameters p,a,T,v^ since V acts by flipping signs while 
taking dual parameters amounts to interchanging p and v) we have a similar result 
for 0(7; A)A(7; A). Combining these results and using Hfi.l()|l now yields the desired 
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symmetry of S for the restricted parameter set. These extra conditions on the 
parameters can be removed by analytic continuation (compare with the proof of 
Theorem ESI). □ 

Remark 6.8. The symmetries described in Lemma |6.7I should be compared to the 
D4 symmetry of R. Note that for R only an 6*3 C W4 symmetry holds 

trivially from its integral representation H3.2|l . where S3 acts by permuting 71, 72, 
and 73. 

Let us now consider the asymptotics of S, compare with the asymptotics 
of R. 

Lemma 6.9. Let wi, ^2 G K>07 7 G C'*, and A G C \ M such that C > 0, where ( 
is given by H5.5|l . Then 

Sir, X, A) = o(e"(|s^WI-5R(7o)— 

for 3?(x) — *■ ±00, uniformly for 3(a;) in compact subsets of (— CiC); where a — 
2Tr/wiW2- 

Proof. Under the parameter restrictions as stated in the lemma, S does not have 
x-independent poles (see Lemma [6.3|l and the integral representation H5.7|l for 'ip 
holds. 

In view of H6.10|l and the asymptotics 

(6.11) , / , = ©(e^F"^"^) 

A(7;a;) 

for 5ft(x) — )■ ±00, uniformly for 3(a;) in compacta, it suffices to prove 

(6.12) V'(7;a;,A) = C'(e"(l»(^)l-"')I3^(^)I) 

for 5R(x) — > ±c», uniformly for 3(a:) in compacta of (— C). The asymptotic for- 
mula H6.11|l follows directly from the estimates H2.10|l and (|2.11|l for the hyperbolic 
gamma function. 

Note that it suffices to prove (16.12(1 for ^{x) — > 00 since 

(6.13) V(7;a;, A) = ■(/'(7;-a;,-A) 

where 7 = {a,p,v,T) (in the 7^ notation, 7 = (—70,71,721^73))- Equation (|6.13|) 
follows by the change of integration variable z ^ — z in H5.7|) and a subsequent 
contour shift. 

To prove (|6.12() for ^{x) 00 we consider the integral representation 1(5. 7(1 of 
tp. We define 

(6.14) e = max{wuW2) + i|5R(A)| + max(|3(p)|, |5(a)|, |3(t)|, \Q{v)\) 

and we consider the division of R in five intervals 

= (-00, -3fi(x) - e), I2 = {-n{x) - e, -di{x) + e), 

^^■^^^ /3 = (-5R(a;) + £,-£), /4-(-e,e), h^{e,oo), 

for 3?(x) > 2e. We write the integral (|5.7|) defining ip as the sum of five integrals over 
Ij (j = 1, 2, . . . , 5) and we bound the integral over each Ij seperately. The intervals 
are chosen in such a way that the estimates ()2.10|) and (|2.11|) for the hyperbolic 
gamma function can be used to bound the integrand over the intervals /i, I3, and 
I5. To estimate the integrals over the remaining intervals I2 and li we use the fact 
that their lenghts are finite and independent of ^{x). For each interval Ij we show 
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that the integral over Ij is e)(e"(|3(^)l-'")l^(^)l) as 3?(a;) oo, uniformly for 

in compact subsets of (— C, C)- As a consequence ip is also of this order. Details are 

given in Appendix 151 □ 

7. Reduction to Askey- Wilson polynomials 

Using an indirect method, Ruijsenaars ITIl Theorem 3.2] proved that R reduces 
to the Askey- Wilson polynomials 1 when the spectral parameter is specialized to 
certain specific discrete values. We now show by a direct calculation that S (I6.1()|l 
reduces to the Askey- Wilson polynomials for the same discrete spectral values. 

Let us first introduce some standard notations for basic hypergeometric series, 
see [5]. For g G C we write 

n-l 

{a;q)n = Y[{l-aq''), 

k 

(ai, a2, . . . , flfc; g)„ = ]^(aj; q)„. 
The q-hypergeometric series is defined by 



s+lVs 



ai, . . . , as+i 



(oi, . . . , as+i;q)k 
^ (bi,...,bs,q;q)k 



provided that either \q\ < 1 or that the series terminates. The Askey-Wilson 
polynomials P are defined as 



rnix;a,b,c,d \ q) = 4^3 



u ^ '^q^q 

ao, ac, ad 



Notice that the g~" term in the above expression causes the series to terminate. 
This implies that r„ is a polynomial of degree n in cosh(27ra;/w2)- Finally, if we 
use the parameter correspondence 

(7.1) a = -e2"T"/™^g, 6 = -e2"^i/™^g, c = -e^"^"/'"^?, d = -e^"^^/'"^?, 
and if we define 

(7.2) A„ = + «7o -I- mwi , 

then the Askey-Wilson polynomials satisfy the Askey-Wilson second order differ- 
ence equation 

C^rn{x;a,b,c,d \ q^) = u(A„; wi, u;2, 7)''n(a;; a, 6, c, d | q^). 

Here we use the Askey-Wilson operator 13.51) and eigenvalue v H3.6|l . 
Ruijsenaars has shown in lllj by an indirect method that 

(7.3) W2,7;x, A„) = rn{x;a,b,c,d \ q^) 

for n £ Z>o, under the parameter correspondence H7.1|l . Similarly we have 
Theorem 7.1. Under the parameter correspondence (|7.1|) we have 
S'(wi,W2,7;a;, A„) = rn{x;a,b,c,d \ q^) 

for n G Z>o. 
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Proof. Without loss of generality we assume that the parameters wi,W2,j,x are 
generic. 

For generic A we can express ^p as an integral 

(7.4) ^h;x,x) = / i{r,x,x,z)dz 

Je 

with I{z) ~ 7(7; a;, A, z) given by 

G(z + a; + A/2 - to/2 ± ir)G(z - A/2 - ™/2 ± ia) 



I{z) 



G{z + x- A/2 + iw/2 ± ip)G{z + A/2 + iw/2 ± iv) 



and with contour C a deformation of R seperating the upward pole sequences of 
/ from the downward pole sequences of /. When A A„, the pole Zk '■= A/2 — 
iw/2 — i(7 — ikwi from a downward pole sequence of / will collide with the pole 
—A/2 + iw/2 — iv + i{n — k)wi from an upward pole sequence of / for < fc < n. 
In order to compute the limit A ^ A„ in H7.4|l . we therefore first shift the contour 
C over the poles at (0 < fc < n) while picking up poles. In the resulting integral 
the colliding poles are on the same side of the integration contour, hence the limit 
A ^ A„ can be taken. 

To calculate the residues of / at Zk we first remark that k consecutive applications 
of the difference equation (|2.3|) yield 

G{z) _ ^k7rz/w2q-k^/2^_^-2-,vz/tU2q.q2-^^^ 



G{z — ikwi) 
Using this equation we can write 

G{z + ikwi + x + A/2 - iw/2 ± zt)G'(z + ikwi - A/2 - iw/2 ± ia) 



I{z) 



G{z + ikwi + x- A/2 + iw/2 ± ip)G{z + ikwi + A/2 + iw/2 ± iv) 

„2k 



X q- 



I ~^{z+ikwi+x-\/2+iw/2±ip) ~ ^ {z+ikwi+\/2+iw I2±iv) 2\ 

. [—e ' ' ' 'q,~e ™2 ^ ' ' q',q )k 



I -■^(z+ikwi+x+X/2-iw /2±iT) -^(z+ikwi-\/2-iw /2±ia) tn 



Using the fact that the residue of the hyperbolic gamma function at z = —iw equals 
()2.9(l . we obtain that the residue Resk of / at Zk equals 

i^wiW2 G{x + X — iw — ia ± iT)G{—iw — 2i(j) 

nesk - 



2tt G{x — ia ± ip)G{X — id ± iv 

2 k 



X q 



. (-e "2^ "q,-e 'q;q )k 



/ — (x-\-X — i(T:tiT) 9 — ^2i(7 9 9 9\ 

(e ""2^ 'q^,e^2 q^,q^;q^)k 

Now we can rewrite S as 

S{r,x,\)^— [^2TiiY,Resk+ f I{z)dz], 

^/wiW2A{-f;x)A{r,X) V ^ Jc ) 

where the contour C is chosen in such a way that all upward pole sequences and 
the poles Zk {0 < k < n) are above C, while all poles in downward pole sequences 
except Zk {0 < k < n) are below C . In this expression the integral J^, I{z)dz has an 
analytic extension to A = A„. Furthermore £'(7; a;, A) is analytic at A = A„, while 
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A(7; A) and Resk {0 < k < n) have simple poles at A = A„. Hence we obtain 

5' 7; X, A„) = hm > Resk 

A^A„ ^wiW2A(7;x)A(7; A) ^ 

^iJP(a; — iii>-i7o) H iH )n 

— C 2 



X 4(/>3 



(e-^(^T0 + '72/3)^-2n.g2)^ 

„-2n „--|5-(2:-™-i7o/l) „- li? (72+73) ^-2n 
9762 ,62 q . „2 „2 



-(2:-iii>+i72/3)„-2n „-;^(7o+7l) „2 

g , e 2 g 



where the notation ^q/i (respectively 72/3) means that there are two terms, one 
with 7o and another with 71 (respectively, 72 and 73). Inserting the parameter 
correspondence (|7.1I) we obtain 

( -2TTx/w2„-l„-2n+2 -27rx/«>2^-l„-2n+2. 2^ 

S('yr A ^ - p27mx/t02„-nif g ° g ^" 

-^i ^iW — " ('„-l„-l„-2n+2 „-lJ-l„-2n+2.„2\ 



X 403 



^-2n^g-27ra7/»2a^g-2W«'25^C-l^-l^-2«+2 ^ ^ ^ 
g-2ira;/iU2|^-lg-2n+2 g-27ra:/t02^-lg-2n+2 '9 '9 



Using Sears' transformation j^i (III. 15)] of a terminating balanced 4^3 series with 
parameters specialized to a = ae~2'^^/™2^ 5 = j^^-2nx/w2^ ^ _ ^-igj-ij^-2n+2^ ^ _ 
ab, e = e-2'^^/"'2c-i5-2»+2^ and / = e-27ra;/i02(;-i^-2n+2 y-^^jg ^j^^ desired 
result. □ 

8. Equality to Ruijsenaars' hypergeometric function 

We have already seen in previous sections that Ruijsenaars' hypergeometric func- 
tion R and the renormalized formal matrix coefficient 5* have several properties 
in common. They satisfy the same Askey- Wilson second order difference equa- 
tions, they have the same duality property, they specialize in the same way to 
the Askey- Wilson polynomials and their possible pole locations coincide. These 
common properties suffice to show that R and S are equal. 

Theorem 8.1. We have 

(8.1) R{wi,W2,i:,x,\) = S{wi,W2,i;x,X). 

This theorem is equivalent to the following identity between hyperbolic integrals. 
Corollary 8.2. For wi,W2,'3i{'jj) > and \x\, |A|, < w/6 we have 

G{z + x + A/2 - iw/2 ± t(73 - 72)/2)C;(z - A/2 - iw/2 ± ij-fp + 7i)/2) 
G{z + x- A/2 -I- iw/2 ± i(7o - 7i)/2)G(z -I- A/2 + iw/2 ± i{j2 + 73)/2) ^ 
_ G{x + ij2)G{x + i"f3)G{X + ^72)G(A + i^a) 
~ G(x + iio)G{x - iji)G{X + i%)G{X - i^i) 

G{z±x + no)G{z ± A + i%) 

3 

c G{z + iw) Ylj^i G{z + ijQ + ijj + iw) 

where the contour C is the real line with a downward indentation at the origin. 

Proof. The proof consists of inserting the integral representations of R and S in 
H8.1(l . See H;-{.2|l for the integral representation of R, and (|5.7|l , l|6.10() for the integral 
representation of 5. □ 
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In order to prove Theorem 18.11 we first consider the Casorati-determinant of S 
and R in the iwi direction. 

Lemma 8.3. The Casorati-determinant 

5(7; z + iwi/2, X)R{j; z — iwi/2, A) — S{'j; z — iwi/2, X)R{'j; z + iwi/2, A) 
of S and R in the iwi direction is identically zero. 

Proof. We suppress the A and 7 dependence of 6{z) whenever this does not cause 
confusion. We prove the lemma for generic parameters wi,W2 & ^>0j 7 G R*, and 
A e J7 \ M, under the condition W2 > 2^ + 2|SJ(A)| + 3wi, where U is an open subset 
such that the asymptotics (|3.11|l of R hold for A G C/. 

A simple calculation involving the Askey- Wilson difference equations satisfied 
by R and S (see (|3.7|l and Theorem 16 . 61 respectively^ shows that 

S{z + iw,/2) = '^(^ - ™i/2)' 

where A is defined by (|3.4|) . Since the function 



T{z) = sinh(27rz/w2) 



G{z — ijj — iwi/2) 



^^^Giz + ijj+iwi/2) 
satisfies the same difference equation, we conclude that 

. ^ S{Z) 

T{z) 

is an iwi-periodic function. 

We now show that m{z) is an entire function in z. Let us look at the possible 
poles of the Casorati-determinant S(z). By Lemma [6.91 the possible poles of S are 
located at 

±(A+-i7,), (j=0,l,2,3). 

From H3.3|l the possible poles of R are located at the same points. Hence S{z) can 
only have poles at 

±{A+-ij,)±iwi/2, (j = 0,l,2,3) 

Here all sign combinations are possible. Furthermore, using the pole and zero 
locations (|2.8() of the hyperbolic gamma function, we can easily see that the possible 
zeros of T{z) are located at 

±(A+ + ijj + iwi/2), riw2, {j = 0, 1, 2, 3; re Z). 

By the assumption that the parameters are generic, we conclude that m has no 
pole sequences of the form p + ikwi (fc S Z). By the zwi-periodicity of m it now 
follows that m cannot have any poles. 
In the limit '3i(z) — > 00 we have 

T(z) 

uniformly for ^(z) in compacta, in view of the estimates (|2.HJ|I and H2.11|l for the 
hyperbolic gamma function. Here a — 2'k jwxWi as before. 
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Furthermore, using the asymptotics for S (see Lemma l^l^ and for R (see 
we have for 5ft(z) — > oo 

S{z) = C'(e2"(l^(^)l+l'^ol"«')l'^(^)l) 

uniformly for 3(z) in compact subsets of (— C + wi/2, ( — wi/2). Observe that the 
interval {—Q-\-wi/2, ^ — u'i/2) is nonempty due to the conditions on the parameters. 
Combining these two asymptotic estimates we obtain 

(8.2) m{z) = ^ = 0(ga(2|9(A)|-^.)|5R(.)|) ^ q 

T{z) 

for 5R(z) — !■ cx), uniformly for ^(z) in compacta of (— C + u'i/2, — wi/2). 

The asymptotics of m(z) for 5R(z) — ^ —oo can be obtained in a similar way and is 
also given by (|8.2|l . Combining the asymptotics with the fact that m(z) is analytic 
and zwi-periodic we conclude that m{z) is bounded on C since — wi/2 > u'i/2. 

For these parameters we conclude by Liouville's theorem that m(z) is constant. 
In fact, by the asymptotic expansion H8.2|l . m is identically zero. We can now 
extend this result to all values of the parameters by analytic continuation, which 
proves the lemma. □ 

Proof of Theorem \8.1[ Consider the quotient 

Q{r,x,\) = — --. 

By Lemma |8. 31 Q is an iwi-periodic mcromorphic function in x. Since Q is sym- 
metric in wi and W2 (for both R and S are invariant under interchanging wi and 
W2), Q is also iw2-periodic. If we choose wi,W2 > such that wi/w2 ^ Q, then 
the set {kwi + lw2 | /c, ^ G Z} is dense on the real line, hence Q(7; x, A) is constant 
as meromorphic function in x. Analytic continuation (in wi, W2, and 7) allows us 
to extend this result to all possible values of wi and W2 in C+ and 7 G C*. 

By the duality properties of R and S (see (|3.8ll and Theorem 16 . 51 respectivelv) . 
we have 

(5(7; a;. A) = g(7;A,a;). 

This implies that Q is also constant as function in A. 
In particular we have 

^^ S{wi,W2,j;x,\q) 
R{'wi,W2,j;x,Xq) 

with Ao given by (j7.2|l . By Theorem 17.11 we have 5'(wi, W2, 7; a;, Ao) = 1, and by 
(17. 3|) we have i?(wi, W2, 7; x, Aq) = 1. Hence Q = 1, as desired. □ 



Appendix A. Eigenfunction of TT\{Yp) 

In this appendix we give the explicit calculation to rewrite the eigenvalue equa- 
tion TT\{Yp)f — HT{p)f as the first order difference equation H5.3|l . 

Using the explicit expression 1)5. l|l of Yp, the eigenvalue equation becomes 

^q-'/\xiE)f + ^q-'/\,iFK)f ^^{n,{K 1))/ = ^^^^f. 
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By the explicit definition fLcmma l4.5|l of tta we obtain 



q-q 

i 



q-q- 



-(/(z + zu;i)-/(z))=^^^/(z). 



Multiplying hy q — q~^ and rearranging the terms yields 

which is equivalent to 

^- -J cosh(7rj/2 + 27rz/w2 — 7rit(;i/(2ii;2) + 7rA/ty2) + cosh(27rzT/w2) ^.^ 
cosh(— 7ri/2 + luzjuji + ixiw\j (^w-i) — ■k\/w2) + cosh(27rip/i(j2) 
Replacing the variable zhy z — iwi/2 we can now rewrite the latter equation as 
f(z + iwi/2) 

cosh(7ri/2 + ^ixzjwi — Z-Kiwi/ {2w2) + 7tX/w2) + cosh(27riT/u'2) 



cosh(— 7ri/2 + 27rz/w2 — Triwi / {2w2) — ttX/'W2) + cosh(27rip/w2) 
cosh(^(z + A/2 - 3iwi/4 + ^2/4 ± ir)) 



f{z~iw,/2) 



cosh(;;^(z - A/2 - iwi/4 - itU2/4 ± ip)) 
'z + X/2-3iw/2±iT)) 

f{z-iwi/2), 



fiz - iwi/2) 



cosh(;^(z + A/2 - 3iw/2 ± ir)) 



coshi^iz~X/2-iw/2 + ip)) 
where we used the j7r-antiperiodicity of the hyperbolic cosine in the last equality. 

Appendix B. The limit behaviour of V' 

In this appendix we give the details on the calculation of the limit behaviour of "0, 
cf. the proof of Lemma 16.91 Throughout this section we assume that wi,W2 G R>0) 
A ^ R, 7 e C^, and that C > (with C given by (ESJ). We prove that 

(B.l) iPir, X, A) = C)(e"(|a(A)|-»)|3i(:^)l) 

for 5ft(x) — > 00, uniformly for 3(a;) in compacta of (— C)C)- -^s explained in the 
proof of Lemma 16.91 we prove (|B.1() by splitting R in five intervals and bounding 
the integral representation (|5.7|) of over each interval. 

B.l. Preparations. Let us first define a function K by 

G{z + X/2~^w/2±^a) 

^' ^) ^ w^xj2TJ^;;j2±^^) ■ 

The integral representation H5.7|) for ^ can then be written as 

(B.2) ip{'-f;x,X)= / K{z + X, X,T, p)K{z, —X,a,v)dz. 

Jr 

The behaviour of K in the limit z — > ±00 is controlled by 

K±{z,X,a,b) = e^M<^-^^)-a-/2+b'm^ 
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Explicitely, for fixed a, 6, and A we have 

(B.3) K{z, A, a, h) = X±(z, A, a, h)eS^'^^-"-''^ 

for ±5R(z) > max(wi,W2) + |3?(A)|/2 + max(|3(a)|, where 

(B.4) |g(z, A, a, 6)1 < c'(!3(2))e-"'"'"(«'i'"'^)l5^(^)l/2, 

with C depending continuously on ^(z), cf. H2.10|l and (|2.1HI . 

B.2. General estimation scheme. Let e and the intervals Ij {j G {1, . . . , 5}) be 
defined as in (|6.14|) and Ht).15|) . We only consider the asymptotics for ^{x) oo. 
Assume that > 2e, causing the intervals to form a partition of the real line. 

We write the integral (|B.2p defining -0 as 



5 

(B.5) V(x)=^^,(x), 

i=i 

where 

'ipj{x) = / K{z + X, X,T, p)K{z, —X,a,v)dz 
Ji, 

for j e {1,2,..., 5}. We bound these intergrals using l|B.3|l (if one of them is 
applicable for the interval at hand). 
For J = 1 we have 

/.-!R(x)-e 

Mx) = / -ft:_(z + x,A,r,p)ii:_(z,-A,a,w)e9i(^+^'^)dz 

J —OO 

iaXx —awx ia{p'^ -^v'^ —t'^ —a'^) /2 I 2awz-\-gi{z-\-iQ{x) ,x) j 

O t> O I O lAiA/ 

J — oo 

= 0(e-"('a(^)+«')SR(^)) 

for 3?(.t) oo, uniformly for 3(.t) in compacta of (—CC)- Here gi{z,x) — 
g{z, A, T, p) + g{z — x, —A, a, v) which satisfies an equation like ljB.4|) for z < — e 

\gi{z + i'^ix),x)\ < Ce-""""("'i^"'^)l^(^)l/2. 

where the constant C is independent of 3(a;), because 5(x) is bounded. In particu- 
lar, gi(z+iQ{x), x) is uniformly bounded for z € (— oo, — e) and a; G {z G C | 3fJ(z) > 
2e,|3(z)|<C}. 

Likewise we have for j — 5, 

/oo 
K^{z + X, A, r, p)K^{z, -A, cr, u)e3=('''^^dz 

/oo 
g-2aiLiz+g5(2,2;)|^2: 

^ 0(-ga(3(A)-«,)!ff(a;)>) 

for 3?(x) — > oo, uniformly for 9(a;) in compacta of (— C,C)- Here 55 is a function 
which satisfies a bound like HB.4|I for z > e, cf. the previous paragraph. 
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For j = 3 we need to be a bit more careful. First observe that 

where (73 = 17(2: + x, A, r, p) + g{z, —A, a, v) is bounded on z € (— 5f (a;) + e, — e) by 
^g-ami„(»i,»2)min(-2:,z-Si(a:))/2^ j^^j j^g^^^g ^j^g Constant C itsclf. Therefore we 

have 

^-SR(2^)+e 

for 5ft(a;) — > 00, uniformly for Q{x) in compacta of (— CiO- Here we get the final 
approximation by evaluating the integral and using that 9(A) ^ 0. 

For j = 4 we cannot use (|B.3p for the entire integrand. However we still have 

^4(x) = J K+{z + X, A, r, p)K{z, -A, a, i;)eS'''^'"^)dz 

= e'"^(""~^) j K+{z,X,T,p)K{z,-X,a,v)e3^^'''''Uz 

^ 0(ga(a(A)-t«)S(a;)^ 

for 3fJ(x) — > 00, uniformly for in compacta of (— Here 54 is a function 

satisfying the bound 54(^,0;) < Ce-""""("'i'"'^)('^(=")-'^)/2 < C, for z e [-e,e] and 
3fi(a;) > 2e. 

Finally for = 2 we have in a similar way 

V'2(a;) = / i^(z + i3(a;),A,T,p)/^_(z-3fJ(x),-A,a,u)e»^(^'^)dz 



^ g-^aSff(x)(-A-»«,) y ii:(z + 5(a;),A,r,p)ii:_(z,-A,CT,u)ef^(^^^)(iz 

for 3fi(x) — > 00, uniformly for in compacta of (— CiC): where 52 is a bounded 
function, cf. the previous paragraph. 

By (|B.5(I we conclude that the asymptotics (|B.ip for holds, as desired. 
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